RBRC-562 
UMDOE40762-351 

Threshold Resummation for Higgs Production in Effective Field Theory 
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We present an effective field theory approach to resum the large double logarithms originated 
from soft-gluon radiations at small final-state hadron invariant masses in Higgs and vector boson 
(-y*,W,Z) production at hadron colliders. The approach is conceptually simple, independent of 
details of an effective field theory formulation, and valid to all orders in sub-leading logarithms. 
As an example, we show the result of summing the next-to-next-to-next-to leading logarithms is 
identical to that of the standard pQCD factorization method. 
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1. Introduction. In hadron colliders such as Tevatron 
and LHC, the rates on Higgs boson and Drell-Yan pair 
production demand reliable perturbative quantum chro- 
modynamics (pQCD) calculations. When the final-state 
invariant mass of hadrons is small, a fixed-order pQCD 
calculation yields large threshold double logarithms in 
the coefficient functions 



ln TO_1 (l-z) 



(1-z) 



(m < 2k) 



(1) 



which must be resummed to all orders in a s , where 1 — z 
is the fraction of center-of-mass energy of the initial par- 
tons going into soft radiations. In moment space, these 
large logarithms appear in the form, ln m N, where 
N = N exp(7#) with N, the order of moment and ^e, the 
Euler constant. In the past decade, a standard method 
based on pQCD factorization has been established to per- 
form the resummation Q, 0, Q , and it has been carried 
out to next-to-next-to-next to leading l ogar ithms (N 3 LL) 
for the above processes 0, II H 0, B 1, [l3 • 

In this paper, we present an alternative, effective field 
theory (EFT) resummation of these large threshold log- 
arithms. The EFT approach is conceptually simple, and 
is readily extended to other processes. It is motivated 
by the recent development of soft-collinear effective field 
theoryJTljll and its applications to threshold resumma- 
tion |l3l The major steps described here are stan- 
dard in EFT methodology and are similar to those in 
Refs. 0, 0| . However, there are also significant dif- 
ferences: First, we work in the limit of Q — Mh (Higgs 
mass) going to infinity first (twist expansion) , and z — * 1 
subsequently. The kinematic region (1 — z) <~ Aqcd/Q, 
where there might be subtleties in formulating an effec- 
tive theory |15| . is avoided. Second, our approach uses 
the full QCD results in the soft-collinear limit, and thus 
an actual formulation of an effective lagrangian is unnec- 
essary. Finally, our resummation formula is valid to all 
orders in sub-leading logarithms and fully equivalent to 
the standard resummation. As an example, we demon- 
strate that the EFT resummation to N 3 LL agrees com- 



pletely with the latest result in the literature using the 
traditional approach 0, [^j . 

Let us consider the standard model Higgs production 
in hadron- hadron collisions. An almost identical discus- 
sion applies to Drell-Yan production of 7*, W and Z 
bosons. At high energy, the Higgs cross section is domi- 
nated by gluon-gluon fusion through the top quark loop. 
For a large range of Higgs mass, the top quark can be 
considered as heavy and can be integrated out first (l6| . 
The effective lagrangian density for Higgs production is 



C = --C^MufiR) d> CTG^r) , 



(2) 



where 4> is the scalar field, G^ v is the gluon field strength, 
Cff, is the effective coupling |17| , and [Ir is a renormaliza- 
tion scale. We focus on the kinematic limit in which the 
final-state hadron invariant mass (1 — z)Mh is small in 
the sense that (1 — z) <C 1. The scale (1— z)Mr, however, 
is still perturbative and is in principle much larger than 
\/Aqci>Mh- In the above kinematic limit, the pQCD 
process is dominated by soft and collinear gluon radia- 
tions. 

Renormaliztion scale \xr is arbitrary in principle, but 
we will set it to Mh in the following discussion. The 
process then has three independent scales: Mh, (1 — 
z)Mh and (If, last of which is related to factoring the 
collinear divergences into Feynman parton distributions. 
To calculate the Higgs production cross section reliably 
in the threshold region, we shall study the physics at 
different scales separately using EFT techniques [Tflj . 

2. Physics at scale Mh, and between Mh and (1 — 
z)Mh- Momentum scale Mh is confined to the gluon- 
Higgs-gluon vertex region. To account for their contri- 
bution, we introduce the scalar current Jqcd(Mh) = 
G^ v G fn,{M h) and make the operator expansion, 



Jqgu(Mh) = C g (M H /n, a a (fi))Jeg(fi) 



(3) 



where coefficient function C g contains the perturba- 
tive contribution between momentum scale Mh and /x. 



2 



J e g(/i) contains the soft and collinear contributions be- 
low scale n. To calculate the coefficient function, we take 
the matrix element of the above equation between gluon 
states — the left-hand side defines the gluon form factor 

Fg(M H ). 

In dimensional regularization and minimal subtraction 
scheme, the gluon form factor can be calculated as a 
power series in the renormalized strong coupling, 



F g (M H ) = Z G 2(M H )J2(a s (^f 



1=0 



F, 



(0 



(4) 



where Zq2 is a renormalization constant, F g 1 ^ contains in- 
frared 1/e-poles. Following y|, we introduce o s = a s /47r 
as expansion parameter for the form factor calculated at 
higher order. The infrared divergences can be isolated by 
the following factorization, 

F g (M H ) = C g {M H /n, a a (jj)) S g [M„/^ a s (/z), 1/e) 

where S g contains the pole contributions only and can be 
regarded as the gluon matrix element of effective current 
J e ff(/u). Expanding the coefficient function at \i = Mh 

as C g (l, ot s {M H )) = J^i a l{^H)Cg l \ and using the form 
factor calculated to the second order in a s 0| > we obtain 



= 7C A ( 2 
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where £„ is the Riemann zeta- function, nf is the number 
of massless quark flavor, Ca — N c and Cf = (N% — 
l)/2N c and N c = 3 is the number of color. 

The physics between Mh and M#(l — z) can be ac- 
counted for by integrating over the running scale fj, be- 
tween them. This requires an anomalous dimension, 



. . dhiCq dhiS, 
7l!s(a5 ) =M __ = _ M __ 



9 — 



E 



i (i) 

a »Ti,s 



(7) 



A simple way to calculate this is to start with the follow- 
ing representation of the form factor [lflj . 



f( 



K g {a 8 {n),i 



(8) 



+ G g (l,a s (^,e),e) + 



1 d\ 



A g (a s (\fi,e) 



where A g is the cusp anomalous dimension of Wilson 

lines in adjoint representation |2jJ, A g = Y\i a \A$ , and 
has been calculated to three-loops recently |2lj| . The soft 
function K g contains only the infrared pole terms and 



can be constructed from A g through [^dKgjdp? = —A g . 

The perturbative function G g — a \Gg contains both 
finite terms and those vanishing when e — > 0. According 
to Ref. j2^|, G g has a simple representation 



Gf = 2B« - 2i#_i + /« 



AG g l) 



(9) 



for i = 1, 2, 3, where AG g can be constructed from 
those vanishing terms in lower-order G g . -62,9 = 
J2i a s-^2*g i s * ne coefficient of 5(1 — x) term in the 
gluon splitting function, with B^ g — 11CU/3 — 2n//3, 
etc j2l|. The QCD /3-function is defined as (3(a s ) — 
—d\n.a s /d\n.^ 2 = /? a s + P\a 2 s + .., with /3 = 11(7^/3 — 
2n//3. The functions /g 1 " 1 are universal in the sense that 
the corresponding quark expressions are obtained by re- 
placing the overall factor of Ca by Cf- 

The pole terms in Eq. (8), the S g terms, will be used 
to calculate the anomalous dimenstion by employing the 
last two relations of Eq. (7) . After taking into account the 
renormalization of the strong coupling constant a s (/j), we 
get the anomalous dimension of the gluon current as 



7 fi = Af In (M^/ M 2 ) + B« + 2ifl 



where 
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(10) 



(11) 



Since AW, -B^g, and /g 1 " 1 are known to three loops |22l 
|23| , the anomalous dimension is now known to the same 
order. 

Using the evolution equation, one can summarize the 
physics between scales Mh and /ij ~ (1 — z)Mh in 
C g (Mff/fi,a a (fJ,)) with a coefficient C g (l,a s (Mff)) and 
running 



V Mi 



C s (l,a s (M H )) 



/ <"Ni) 1 - y 7i, s — I ,(12) 



where the exponent contains Sudakov double logarithms. 

3. Physics at scale (1 — z)Mh- At this scale, one must 
consider soft-gluon radiations from the initial gluon par- 
tons. In principle, one should formulate a soft-collincar 
effective theory to calculate these contributions, as was 
done in Ref. |13| . However, this is unnecessary in prac- 
tice and the result can simply be obtained from a full 
QCD calculation at the appropriate kinematic limit. 

The real emission diagrams without any internal ra- 
diative corrections contain both soft and collinear diver- 
gences. The infrared contributions from diagrams with 
internal radiative corrections (just the pole terms in di- 
mensional regularization) serve to cancel the soft diver- 
gences. The collinear divergences can be factorized into 



3 



a standard Feynman parton distributions. The finite re- 
mainder is the coefficient function (or matching coeffi- 
cient in the sense that we match a product of effective 
gluon currents onto a product of gluon light-cone distri- 
bution operators) which is what we are interested in. To 
see that the relevant physics happens around the scale 
Hi is to take Mellin transform of the coefficient function, 

which contains logarithms of type ln m Mfj / /j, 2 N . If 
\x is set as /// = Mjj/N, the large logarithms disappear. 

Expanding the matching coefficient M-n = 
T,i c i(l J 'i)-M$, the full pQCD calculation H @ in 
the soft limit yields 



M 



(2) 
N 
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(13) 



at scale fj,i = Mjj/N. Note that the coupling constant 
a s is also evaluated at the intermediate scale fij. 

4. Physics between scales fii and (If, and at \xf- In 
the previous section, QCD factorization produces gluon 
distributions at scale fii = Mh/N. We can bring 
the distributions to an arbitrary scale hf using the 
standard DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi) evolution. This introduces an evolution factor, 



(14) 



where the twist-two anomalous dimension 7^ has the 
following large N behavior: 



I2 



N 



-A g \nN + 2B 2 



9 ' 



(15) 



where A g and B 2g are the same as those in Eqs. (8) 
and (9), respectively. The moment of gluon distributions 
introduces a factor, 



g(fi F ,N)g(nF,N) 



(16) 



To simplify the result, the factorization scale fiF is hence- 
forth chosen to be Mjj- 

5. Resummation to all orders in sub-leading loga- 
rithms. Putting all factors together, the Mjj/s-moment 
of the cross section is (s is the total center-of-mass energy 
squared) Q 

a N = a ■ G N (M H ) ■ g{M H , N)g(M H , N) , (17) 

where o~o is a reference cross section and 

G N {M H ) = \C(a s (M H ))\ 2 e^ MH ^ 

xMNiaMy^''^ (18) 

is a pQCD factor, where I\ = 2 4jf7i,g with ji <g — 
7i j9 — 2i/3j_i, is the anomalous dimension for C — C^ x 



C g , and I2 = 2 f£f H -7^72,9- To capture all large IniV, 
we translate the dependence on the intermediate scale 
a s (/Ar) to a s (M H ), using 



MjvMm/)) = M N {a s (M H )) e h 



(19) 



where I 3 = -2 f ^AB X , and AB 1 is defined as 
ABi = -fi(a s )d\nMN /dhia. s . 
The final form of resummation is 



G N (M H ) = TiasiMH^e 1 ^^ 1 ^ 



(20) 



where T = \C(a s (M H ))\ 2 M(a s (M H )) depends only on 
a s {Mn)- In fact, it is just the standard full-QCD coef- 
ficient function in the soft-gluon approximation, without 
the large logarithms. / = /1 + I 2 + I3 is a function of 
A = (3q In Na s (MH) and a s (Af^) with all leading and 
sub-leading large logarithms resummed. 

The above result can be related to the conventional 
expression if one writes 



/ = / A +hiAC* , 



(21) 



where 



dz- 



1 - z 



(l-z) 2 M 2 H A „i 



Ml 



+D g {a s {{l-z) 2 M 2 H ))} 



(a s (M 2 )) 
(22) 



and AC is just a function of a s (Mtf), serving to cancel 
the non- logarithmic terms in 2a- Using similar methods 
as the ones in Ref. |7j, it is a matter of some technical 
steps to get, 



D g (li 2 ) 



AC = 



2(B hg + ABi + 2B 2 , g ) 
-d a T 2 (d as ) [AA g {a s ) - d a D g {a s )\ 
T 2 (d as )[AA g {a s ) - d a D g {a s )\ , (23) 



where r 2 (e) = 1/ 1 e 2 [l-e~-<^T{l-€)] = -C 2 /2-C 3 e/3+... 
and da 3 = 2(3(a s )a s d/da s . The above equations are 
our main result connecting the EFT resummation to the 
coventional approach, valid to all orders in leading and 
sub-leading logarithms. Similar results have been ob- 
tained for deep-inelastic scattering and Drell-Yan pro- 
cesses |24j . 

6. Results up to N 3 LL and Conclusion. The function 
/ has a perturbative expansion in o. s (Mh)i 



I = hxN 9l {\) + g 2 (X) + a s (M H )g 3 (X) 
+a 2 s (M H )g 4 (X) + ... , 



(24) 



where <?i(A) can be found from Ref. 0,0. They are 
functions of A g and D g : g\{\) sums over the leading 

logarithms, depending on A g , and 52(A) sums over next- 
to-leading logarithms, depending on A g 2 ^ and Dg 1 ', etc. 
The D g ^ coefficients can be solved iteratively from Eq. 
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(25) 



and so on. 

As an example to demonstrate the equivalence to 

(3) 

the conventional formalism, we calculate D g using the 
known result f { g 3) [H, M^<^ ( Eq . (13)), and A g 1] ' (2) 
plf. The answer is, 
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(26) 



which agrees completely with the recent calculations |9j, 



To conclude, we have presented an effective held the- 
ory method to resum large threshold double logarithms 
in standard model Higgs production. The approach is 
simple conceptually and uses the full QCD calculation in 
the soft limit. The result is valid to all orders in leading 
and sub-leading logarithms, and reproduces the known 
answer to N 3 LL order. 

It has been shown that the inclusion of threshold re- 
summation effects helps to reduce the theoretical uncer- 
tainties in the prediction of Higgs production rates at 
hadron colliders 5] . It will be interesting to see if we can 

(3) 

get even better results when the newly calculated D y g ' 
term being included. We note that the other sources of 
theoretical uncertainties, such as those stemming from 
heavy quark loop approximation |16L , and those from 
parton distribution parameterizations, need to be consid- 
ered in a detailed phenomenological studies of the Higgs 
production. 
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